Enhanced observability of quantum post-exponential decay using distant detectors 
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We study the elusive transition from exponential to post-exponential (algebraic) decay of the 
probability density of a quantum particle emitted by an exponentially decaying source, in one 
dimension. The main finding is that the probability density at the transition time, and thus its 
observability, increases with the distance of the detector from the source, up to a critical distance 
beyond which exponential decay is no longer observed. Solvable models provide explicit expressions 
for the dependence of the transition on resonance and observational parameters, facilitating the 
choice of optimal conditions. 
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I. INTRODUCTION 

Exponential decay results when the decay rate de- 
pends linearly on the surviving population. It is ob- 
served in many quantum systems, in virtually all fields 
of physics (particle, nuclear, atomic, molecular, and con- 
densed matter), but its microscopic derivation from the 
Schrodinger equation is not as straighforward as in clas- 
sical physics, due to initial state reconstruction [J, and 
deviations are predicted at both short and long times 
0, H, The Zeno effect, in particular, is associated 
with the short time deviation. The deviations at long 
times are less often discussed, and constitute the central 
topic of this paper. 

On the experimental side, the post-exponential regime, 
which is typically algebraic, has been elusive. There 
have been many attempts to produce evidence of post- 
exponential decay, with little success [^, 0, B ■ It has 
been argued that repetitive measurements on the same 
system, or simply the interaction with the environment 
would lead to persistence of the exponential regime to 
times well beyond those expected in an isolated system 
01 [13, [HI ■ Nevertheless a recent measurement of transi- 
tions from exponential to post-ex pon ential decay of ex- 
cited organic molecules in solution [Tj] appears to contra- 
dict this expectation, and has triggered renewed interest 
in the topic. 

Aside from the non-trivial challenge of understanding 
and extending these experimental results, there are many 
reasons for studying post-exponential decay: Winter, for 
example, argued that hidden-variable theories could pro- 
duce observable effects in samples that have decayed for 
many life-times [l3j ; more recently Krauss and Dent have 
pointed out that late-time decay may have important cos- 
mological implications [ij]; decay at long times is quite 



sensitive to delicate measurement and/or environmen- 
tal effects, so it is a testing ground for theories of these 
processes; at a fundamental level, the form of long time 
deviations from exponential decay might distin guis h be- 
tween standard, Hermitian quantum mechanics 15| , and 
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modifications with a built-in microscopic arrow of time 
[l^ . [l7]; from a practical perspective, Norman pointed 
out that post-exponential decay could set a limit to the 
validity of radioactive dating methods 0; we have also 
argued that the deviation could be used to characterize 
certain cold atom traps Indeed, due to technologi- 
cal advances in lasers, semiconductors, nanoscience, and 
cold atoms, microscopic interactions are now relatively 
easy to manipulate: this makes decay parameters con- 
trollable, and post-exponential decay more accessible to 
experimental scrutiny and/or applications. For example, 
under appropriate conditions it could become the domi- 
nant regime and be used to speed-up decay, implementing 
an Anti-Zcno effect [19]. In addition, recent experiments 
on periodic waveguide arrays provide a classical, elec- 
tric field analog of a quantum system with exponential 
decay pol . l2l[ . These experiments may also reach the 
post-exponential region in a particularly simple way. 

Most of the proposals to enhance post-exponential de- 
cay, making it more visible by advancing its onset time, 
are based on the idea of observing the survival proba- 
bility for a resonance close to threshold ("small Q- value 
decay") [H, [H [H, [H, [13, HI , or more precisely, with 
energy release comparable to resonance width. This is a 
rare configuration in naturally decaying systems, such as 
radioactive isotopes , or it leads to difficult measure- 
ments [13, [H, [1^ , and realizing it in artificial structures, 
although possible in principle, is still a pending task [23 |. 
The only realization so far may be the experiment on 
organic molecules in solution, mentioned above [T^. but 
this is a rather complex system and its analysis from 
first principles has not yet been performed, so the true 
mechanism behind the observed data remains to be con- 
firmed. Another idea is to use the escape of interacting 
cold atoms from a trap, in the strongly interacting Tonks- 
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Girardeau regime, to enhance the signal [13] ■ Again, the 
corresponding experiment has not yet been done. We 
also mention the proposal of Kelkar and coworkers [l^ 
to deduce, rather than observe directly, the long-time de- 
cay characteristics from phase-shift data. 

In the present paper we propose a quite different 
strategy with the important advantage of simple im- 
plementation. We show, using solvable decay models, 
that by increasing the distance of the observation point 
from the source, the transition from exponential to post- 
exponential decay occurs at higher probability densities, 
and thus becomes more easily observable. Figure [T] shows 
the probability density versus time, for three positions, 
and illustrates the density increase at the transition as 
the observation position x is increased. (The details of 
the model will be explained in Sec. II.) The observabil- 
ity of the transition may additionally be enhanced by 
amplifying the signal with a macroscopic Bose-Einstein 
condensate wavefunction as explained in the final dis- 
cussion. The combination of enhancing factors makes it 
possible to predict the observability of the transition for 
longer lifetimes than those implied by the above men- 
tioned "small Q-value" condition. 
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FIG. 1: (Color online) Probability density versus t for a fixed 
value of fco = 1 — 0.5j. Three detector positions x have been 
considered: x = 1.5 (red dotted line), x — 0.4 (blue solid line), 
and a; = 0.1 (green dashed line). The dimensionless variables 
and the model are explained in Sec. II. 

Finally, note that when the distance is too large the ex- 
ponential decay regime disappears, and with it the tran- 
sition, so that only non-exponential decay should be ob- 
served (28| . 

In Section II a simple decay model is introduced and 
solved exactly. Section III provides approximate expres- 
sions useful to define the transition from exponential to 
post-exponential decay discussed in Section IV. The fi- 
nal section provides examples of possible physical realiza- 
tions and compares the results with those of a discretized 
site model. 



II. EXPONENTIALLY DECAYING SOURCE 
MODEL AND EXACT SOLUTION 

There is a long and fruitful tradition of analytically 
solvable models for studying diverse aspects of single 
particle quantum dynamics. Closest to the model pro- 
posed here are several "source models" to examine quan- 
tum transients (propagation of signals, diffraction in 
time, forerunners) and tunneling dynamics of evanes- 
cent waves Many works have been carried out 
within this approximation to study neutron interferom- 
etry Isol. [s^. atom- wave diffraction [s^l, tunneling dy- 
namics [29l l33l [3^. [bsI. [s^ . [stI [ssI . absorbing media [s^, 
different aperture functions [3Cll.l32. 40], and to model the 
dynamics of a guided atom laser |41j] . 

In these "source models" , unlike the usual initial value 
problem, the wave function is specified for all times at 
one point of coordinate space. (The connection between 
initial value problems and source boundary conditions 
has been discussed in [i^l). A simple, solvable ID case 
corresponds to switching-on the source suddenly at T = 

0, with a well defined "carrier frequency" Oq > 0, 

*(o,r) = e(r)e~'"'''^, (i) 

zero external potential, ^ = 0, and a vanishing wave ev- 
erywhere (for all X) before the source is switched on, i.e. 
for T < 0. The solution of the Schrodinger equation for 
X > is identical to that for the "Moshinsky shutter" , 

1. e. a wave created by the sudden opening of an opaque 
shutter which totally reflects an incident plane wave with 
reflection amplitude r = 1 [3^] prior to T = 0. This 
solution shows the "diffraction in time" phenomenon of 
a fairly well-deflned wavefront advancing with velocity 
Koh/m, where 

Ko = y^2mno/h (2) 

and m is the mass of the particle. Negative values of 
correspond to particle injection with energy below the po- 
tential threshold and may also be considered: interesting 
transients are produced in this way and, at long times, 
a stationary evanescent wave which penetrates into the 
classically forbidden region [2^, [33| . Here we shall work 
out a case that has not yet been examined, namely, a 
complex carrier frequency fio = f^oi? + i^oi- 

We assume a negative imaginary part, ^Iqi < 0, so that 
the density at the source point decreases exponentially, 

|1'(0,r)|2 = e~2|Ooi|T^ (3) 

The real part floR is chosen to be positive so as to gen- 
erate travelling waves released from the source, and thus 
mimic the behaviour of resonance decay above threshold, 
in the outer region (outside the interaction region that 
holds the resonance).^ 



There are no resonances below the energy threshold, so we disre- 
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This provides a simple (minimal) solvable model for 
an exponentially decaying system with a long lived reso- 
nance, while avoiding a detailed description of the inter- 
action region where the decaying system is prepared. For 
many applications, and in particular during the exponen- 
tial decay of the initial state, these details are unneces- 
sary in the outer region, which in our case is represented 
by the positive half-line. We have also checked, see the 
final discussion, the broad validity of the results obtained 
using a different, complementary, norm-conserving, dis- 
cretized lattice model. In it, the particle is prepared in 
the first site, which plays the role of the initial trap from 
which it escapes. This discretized model can be realized 
in periodic waveguide arrays, making use of an electric 
field analog of the quantum system [23, HH . 



A. The dimensionless Schrodinger equation 

The dynamics of a particle in free space, V = 0, 
given by the Schrodinger equation 



^5* 



2m dX^ ■ 



(4) 



Let us first reduce the number of variables by introduc- 
ing dimensionless quantities for position, time and wave 
function, in terms of some characteristic length L: 



X 

t 



X/L, 

Th 

^ ^{X,T). 



(5) 
(6) 
(7) 



L may be chosen for convenience and we take it as K^^ , 
where Kq = (2mf2o/^)"^^^ = Kqr + iKqi is the di- 
mensioned complex wave number and Kq^ its real part. 
The branch cut is drawn just below the negative real 
axis so that for fio/ < the imaginary part of K^i is 
negative, Kqi := lm{KQ) < 0. Defining dimensionless 
wavenumber and carrier frequency as ko — KqL, and 
ujQ = ^lo2mL^/h = fcg, this election fixes the real part 
of the dimensionless wave number to be unity, fcoij = 1, 
and the dimensionless dispersion equation becomes 



(8) 



where — 1 < koj < 0. In terms of these new variables, 
the Schrodinger equation ([4]), is now 



dt 



(9) 



gard the case f2o-R < Oi even though it could be treated formally 
in parallel with the Qq-B > case. The evanescent case SIq-R < 
might represent approximately the decay of a resonance in a tun- 
nelling region of space far from the particle release to the outer 
region. The main differences with respect to the travelling case 
will be pointed out in a later footnote. 



and the dimensionless continuity equation is 
dp{x,t) dJ{x,t) 



dt 



+ 



dx 



where 



p{x,t) = |^(x,t)|2. 



J{x,t) = 2Im 



= 0, 



d^{x, t) 
dx 



(10) 

(11) 
(12) 



are the dimensionless probability and current densities. 



B. Explicit solution 

We are interested in waves that decay exponentially in 
time so the Schrodinger equation must satisfy boundary 
conditions at the origin 



V'(o,i) 



B{t), ujoR. > 0, cooi < 0, 



(13) 



and at infinity, where the wave remains bounded as a; — > 
oo. which guarantees that the wave function does not 
diverge at a; = oo. 

To find the solution for a; > we make the ansatz 



(14) 



and determine the function A(lu) from the boundary con- 
dition (fT3| . 



^(0,t) 



, —iujQt 



e{t) 



cLu A(uj)t 



Inverting the Fourier transform we obtain 
Aiuj) 



(15) 



(16) 



The resulting integral for Tjj{x,t) is done most easily in 
the complex k plane, 

k^y^. (17) 

We may rewrite Eq. (jl4p for the wave function as 



V'(a;,t) 



-1 
27ri 



1 



1 



k + ko k ~ ko 



, (18) 

where the branch cut for k = ^/uj (or ko — y/iuo) is cho- 
sen as before, slightly below the negative real axis. The 
integration path runs first from ioo to oo in the complex 
fc-plane, but using Cauchy's theorem it can be deformed 
into the diagonal path r_|_ in the second and fourth quad- 
rants passing above the pole at —ko in the second quad- 
rant, 

^ix,t) = -i / dfce-C^-^-^-) (-^ + , 
^"■^ Jr+ \k + ko k — ko J 

(19) 
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which may be spHt into two integrals of the form 



where 



^-iiak-' +bk) 

3= I dk — - — , a > 0, 



(20) 



The saddle point of the exponent is at {kji = —b/2a, kj — 
0) and the steepest descent path is the straight line kn = 
— {ki + b/2a). These integrals have been studied in [23| 
by contour deformation into the steepest descent path. 
Independently of the pole position at Up = u(kp), the 
result is 



3 = -Z7re''''/4"u;(-up), 



(21) 



where w{z)_:= e ^ erfc(— iz) is the Faddeyeva- (or "w") 
function (43l. |4^. and Up = u{kp), where 



u(fc) 



-(l + z)(fc + 6/2a), 



(22) 



which becomes real along the (straight) steepest descent 
line, and vanishes at the saddle point. 

Now we return to Eq. (I19|) where a = t and b = ~x, 
so the saddle point is at ks = x/{2t): this corresponds 
to the velocity of a classical particle that arrives at x at 
time t having departed from the origin at t = 0. The 
steepest descent path is the straight line kj = — [fc_R ~ 
x/{2t)]. Applying Eq. (pij) . the wave function can finally 
be written as 



ipix,t) 



1 



ikit 



where 



,(+) 



Uiko) = il+l)\j-ko{l~T/t), 



(23) 



(24) 



ui-ko) = -{l + i)\j-ko{l + T/t), (25) 



and the modulus of the complex time 



2fco 



(26) 



generalizes the Biittiker-Landauer traversal time, at least 
formally, in the present context ^29^|. 



Wave function normalization 



N+{t) = / dxp{x,t) 
Jo 



(28) 



is the norm, counting particles emitted to a; > 0. 

As we assume that one particle is emitted as i ^ oo, 
the normalized wave function, denoted by a tilde, is 



'4'{x,t) 



1 



■ip{x,t). 



(29) 



III. ASYMPTOTIC BEHAVIOUR 

The exact solution is useful but not necessarily physi- 
cally illuminating, whereas approximations often provide 
a more transparent physical interpretation. To this end 
we separate each w-function term into the contribution 
of the saddle point k = k^, and a pole contribution. The 
dominant contribution of the saddle is obtained from Eq. 
P9p by setting k = kg in the denominators and integrat- 
ing along the steepest descent path. 



1psix,t) 



1 



1 



1 



2i^ 
[2t 



,(+) 



,(-) 



ikit 



TT (i - l)fco(i2 



(30) 



After the steepest descent path crosses the pole fco when 
Iuiuq^' = 0, its residue ipoix,t) must be added to the 
saddle wave function -05, see Fig. [2l 

„( + )2 



■00 



e = e 



^-iujot^ikox 



(31) 



(Note that the pole at — fco, in the fourth octant, is never 
crossed.) A travelling wave component moving right- 
wards is generated. For fixed t the exponential decay 
is imprinted in e^^"^: because of the negative imagi- 
nary part of fco the wave increases rightwards in coor- 
dinate space until the x value satisfying Imwg^^ = 0. Of 
course the full wave does not show this sharp edge but a 
smoothed one. 

It follows that the total wave function can be approx- 
imated as 

0(x,t) = V,(a;,i) + 0o(a;,i)0[Muo"^^)]. (32) 



For a fair comparison between different fco values, the 
total number of emitted particles must be kept fixed for 
all fco J so the wave fuction must be normalized. The nor- 
malization constant can be obtained from the continuity 
equation (|10|) by integration. 



N+{t = oo) = / dtJ{0,t), 
Jo 



(27) 



To see the conditions for validity of this expression, note 
that the first and the exponential terms in the asymptotic 
series expansion for u;(2;) at large \z\, 



w{z) 



2e-' 



1 X 3 X ... X (2m- 1) 



m— 1 



(2z2 



e[-Im(z)], 



(33) 
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FIG. 2: fc-plane with contour deformations along the steepest 
descent paths (SDP) (the straight lines) at times t2 > ti, and 
saddle points ksi, ks2- At ti, the fco pole has not yet been 
crossed, while at t2 the pole has been crossed and encircled by 
the keyhole contour. The dashed line is the SDP for t ~ oo, 
showing that the pole at —fco is never crossed. 



do reproduce Eq. ([5^ . In our case, large z means large 



,(±) 



Uq ' . Their moduli can be written as 



l„<«|.^^_L±_L + M, (34) 

which take large values in certain conditions, in particu- 
lar at times short and long compared to |t|. 



The transition is generally not sharply defined in the ex- 
act density, and is characterized by interference fringes 
due to the superposition of saddle and resonance contri- 
butions. 

In the following we shall investigate the behaviour of 
the probability density at the transition point tp, when 
the parameters x and/or the imaginary part feg/ of /cq 
are varied. First we fix feg and observe the density as 
a function of t for different x-values. A global 3D plot 
is shown in Fig. [31 in which we chose fco/ to highlight 
the dominance of exponential decay at small x and post- 
exponential at large x. The pole (straight lines) and sad- 
dle densities for three values of x are shown in Fig. HI 
in which x increases from the dotted to solid to dashed 
lines. 
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IV. TRANSITION FROM EXPONENTIAL TO 
POST-EXPONENTIAL DECAY 

The criterion which sets the time scale of the transi- 
tion, is that the moduli of pole and saddle contributions 
to the total density are equal. We thus identify the tran- 
sition as the space-time point where their ratio R{x, t) is 
one. 



R{x,t) 



7r|fcoP[t4 + |r|4-2t2Re(T2 




(35) 
(36) 



(37) 



In addition, the condition Im(uo"^^) > (the pole has 
been crossed by the steepest descent path) must be ver- 
ified. 

Solving R{x,t) — 1 with the assumption t >> |t|, we 
obtain an expression for the transition time tp for fixed 

X, 

„„-Im((^otp-fcox) 

, 12 ■ (38) 



t. 



3/2 _ 



FIG. 3: Overview of the Log probability density in the t, x 
plane, for fixed fco = 1 — 0.15i. 
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FIG. 4: (Color online) Probability density versus time, for 
three x values, at fixed fco = 1 — 0.3j. The crossing point 
between the pole (straight line) and saddle terms defines the 
transition time tp. Dotted (red) line: x = 0.1, crossing at 
tp = 12; solid (blue) line: x = 2.5, crossing at tp = 7: and 
dashed (green) line: x — 12, crossing at tp = 9. 
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(^o/=-0.3) 
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FIG. 5: Probability density versus tp{x) at the transition 
point for fixed ko = 1 — O.Si. The drawing covers x run- 
ning from 0.0001 to 13. The dots correspond to the x values 
considered in Fig. 2] 



Several interesting conclusions can be inferred from 
these figures and the corresponding equations p6l37l38p . 
An important one is that the probability density at the 
transition point increases, as the observation point is 
moved away from the source, improving the possibility 
of observing the transition. A basic reason for this is 
the asymptotic growth of |?Asp ~ x^, whereas the pole 
term density (straight line in Fig. 2]) is basically shifted 
by 5x/2kQji for a shift 5x in the observation coordinate:^ 
in other words, the exponential behaviour is delayed by 
increasing x.^ 

Note also the behaviour of the transition time tp as x 
increases. Initially tp occurs earlier as x increases, but 
there is an x value, x = l/|fco/|, beyond which the ten- 
dency changes, and tp then increases with x until the 
transition is no longer observable, as the onset of the 
exponential decay contribution is below the saddle con- 
tribution. Figure [5] shows the density at the transition 
versus tp. The curve depends parametrically on x, which 
increases from the bottom right corner upwards. The last 
point marks the largest value of x for which the transi- 
tion may be defined by the R — 1 criterion, and such 
that Im(ug~'^'') > 0. This value of x corresponds very 
nearly to the case depicted by dashed lines in Fig. 21 
This maximum x increases with the resonance lifetime 
To = l/|4fco/|, and it is plotted versus |fco/| in Fig. [S] 
The (normalized, see Sec. Ill C[) density at that critical 
point is shown in Fig. [71 We can infer from this figure 
that there is a maximum of the transition point density 
slightly to the right of koi ~ —0.5. This is consistent 



with the condition described by Jittoh [22*1 for purely 
non-exponential decay, 

r/e=l/[fcoVo]=4|fco/| >2, (39) 

where F and e are the (dimensionless) decay rate and real 
energy of the resonance, and we have used the fact that 
the real part koR = 1 in our units. 

We are now better prepared to come back to Fig. [3] 
and appreciate several of its features: note the exponen- 
tial decay at the x = Q edge, the linear delay of the ap- 
proximate wave front when x is increased, the increase of 
the algebraic long-time tail with x, and the characteristic 
transition fringes at intermediate x values, which disap- 
pear for larger x values, characterized by purely non- 
exponential behavior. 
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X 

^ 3.5- 
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FIG. 6: Largest x value where an exponential to post- 
exponential transition is found, satisfying the conditions R = 
1 and Imu}^^ > 0, versus fcoz (normalized waves). 
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^ The onset of the exponential term is also delayed as determined 
by Imn^+' = 0. 

^ This is one of the main differences from the evanescent case 
< 0, for which the pole crossed by the steepest descent 
path is at — fco, and the exponential term is advanced in time by 
increasing x. 



FIG. 7: (Color online) Density of normalized waves, at the 
largest x value where the transition is defined, as function of 
koi. Solid (green) line: density of the approximate solution; 
Dotted (black) line: value from the exact wavefunction. 
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V. DISCUSSION 

The transition from exponential to non-exponential de- 
cay for an isolated decaying system has never been ob- 
served despite the well-established theoretical prediction 
that it must eventually occur. Apart from the difficulty 
of avoiding interactions that may extend the exponential 
region, a fundamental reason is the small value of the 
survival probability at the predicted transition. We have 
shown that observing the density of decay products at an 
optimal distance from the source, significantly improves 
the observability of this elusive transition. Almost 50 
years ago R. G. Newton [2^ gave a classical-trajectory 
argument pointing out that the transition would be af- 
fected by the observation distance. His predicted be- 
haviour and dependence on parameters do not, however, 
coincide with the present quantum results. We have put 
forward a simple decay model that provides analytical 
expressions, facilitates predictions, and allows analysis of 
a proposed experiment through its dependence on res- 
onance and observation parameters (release energy and 
lifetime, distance and time of density measurement). Ex- 
periments may be undertaken using quantum systems 
such as a leaking quantum dot, or escape of ultracold 
atoms from a trap. In both cases, trap and escape pa- 
rameters are controllable, and it is possible to set a time 
zero for preparation of the decaying system. An effective 
ID waveguide (a quantum wire for electrons, or a detuned 
laser for atoms) may channel the escaping particle away 
from the trap. Localized, or spatially resolved particle 
detection can be performed, e.g., with fluorescence or ab- 
sorption techniques. Experimentally accessible low tem- 
peratures and conditions determine decoherence lengths 
and times, large enough to observe the effects described. 
Atom-atom interactions may be minimized by using Fes- 
hbach resonances. As an example, let us consider an ®^Rb 
quasi-one dimensional Bose-Einstein condensate sample 
of 10^ atoms, released to a quasi-lD waveguide from a 
controllable trap that will determine the resonance en- 
ergy and lifetime [45]. Then a commercial CCD camera 
with ~ 3 /im spatial resolution and shot noise limited 
detection with at least ^ 10 atoms/pixel may be used 
to measure the atomic density at the waveguide by reso- 
nant absorption imaging [46j . According to our model a 
sample of 10® atoms with lifetime ~ 400 fis, and release 
velocity ~ 1 cm/s would provide a 9200 atoms/pixel tran- 
sition density at tp ~ 10 ms and at a distance of ~ 100 
/im. These numbers make the observation of the tran- 
sition very plausible. Note that even 10^ initial atoms 
would be close to the shot-noise observation limit at that 
space-time point. An analysis of the effect of atom-atom 
interactions in the condensate is left for a separate study 
but, it is expected that such interactions will not be sig- 
nificant in the transition densities. 

A different physical realization may be possible 
through an array of long parallel waveguides [isj . 



The variation with z, (the longitudinal distance along 
the guide) of the electric field in the n'th guide, plays the 
role of the variation with time of the site amplitudes of 
a quantum tight-binding model f20l . [2]| . Recent experi- 
ments with scanning tunnelling microscopy have allowed 
Longhi and collaborators [2l|] to measure the evanes- 
cent fields along the waveguides, and demonstrate the 
classical-quantal analogy for a periodic parallel array. In 
addition [43| the classical analogue of the quantum Zeno 
effect has been verified. 

If the coupling between adjacent sites, except for the 
first two, is described by a universal hopping constant, 
and the site energies are taken all to be equal, the Hamil- 
tonian is, in dimensionless units. 



H = -A 

OO 

- E 



1>(2| + |2)(1 
|n)(n+ 1| + |n-|- l)(n| 



(40) 



where we take < A < 1 to provide isolation of site 1 
from the rest of the array. 

When only site 1 is occupied at time zero, the analytic 
solution for the amplitude at site n shows a transition 
from exponential to power-law decay, at the time defined 
by equating the exponential contribution to the asymp- 
totic power-law decay expression, dropping a sinusoidal 
oscillation peculiar to this model. 



^3/2 ^ i^ + ^^i^~^)] n^n+l 

V^(i 



v2^3 



(41) 



where = 1 — A^, and 7 — 2A^/a. For small A, 
I a I ~ 1, making 7 << 1, and giving a long lifetime for 
exponential decay. But so long as \a\ < 1, |q;|""'"^ << 1, 
for large enough n, allowing a solution to be found. (The 
oscillations in the power-law decay cause some impreci- 
sion in the transition time, similar to that discussed in 
relation to Fig. [3] for the model of this paper.) Eq. ([4T|) 
is very similar to Eq. (j38|) . and predicts the same phe- 
nomena. This provides support for the general validity of 
our main conclusion, namely, that the optimal transition 
density increases with increasing distance of the detector 
from the trap. 
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